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We investigate the lattice walk performed by a tagged member of an infinite 'sea' of 
particles filling a d-dimensional lattice, in the presence of a Brownian vacancy. Particle- 
particle exchange is forbidden; the only interaction between them being hard core exclu- 
sion. The tagged particle, differing from the others only by its tag, moves only when it 
exchanges places with the hole. In this sense, it is a lattice walk "driven" by the Brow- 
nian vacancy. The probability distributions for its displacement and for the number of 
steps taken, after n-steps of the vacancy, are derived. Surprisingly, none of them is a 
Gaussian! It is shown that the only nontrivial dimension where the walk is recurrent is 
d = 2. 

1. Introduction 

In solids, long-range diffusion of atoms mediated bv vacancies (the so-called 
vacancy mechanism) is an ubiquitous phenomenonElu'tjalj. An atom can perform 
a displacement only if there is a neighbouring vacancy. For low enough vacancy 
densities the motion performed by each of these vacancies can be considered as a 
simple (Brownian) random walk on the crystalline lattice, independent of the other 
walkers. Let us tag one of the atoms. Due to encounters with the vacancies, this 
atom performs some kind of random walk. One question naturally arises: what 
are the properties of this walk? Specifically, we may ask: What is the probability 
distribution to find this atom at a certain site? How long does it take, on the 
average, to get there, etc? One central result of this study is that the probability 
distributions are not Gaussians, being very differrent from the walk performed by 
the vacancy. Since this walk is "driven" by a Brownian hole, we introduced a name 
for it: "the Brownian vacancy driven walk". 

The precise formulation of this problem is stated in Section 2 together with the 
derivation of the complete probability distribution characterizing the couple tag- 
vacancy. In Section 3 we cover some important relations from the theory of lattice 
walks and walks with taboo sites, so that in Section 4, we can explicitely present 
the probability distributions for the displacement and for the number of steps the 
tagged particle made. It is also shown that the only non-trivial dimension where 
the walk is recurrent is d = 2. 

Before we present our analysis, we remark that, after completion of our work, 
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we were informed (Krapivsky) of the study by Brummclhuis and HilhorstO. These 
authors considered the same problem, but only in two dimensions. Our approach 
differs from theirs, in that we considered systems in arbitrary dimensions: d. By 
leaving d as a parameter in the theory, we see the uniqueness of d — 2 systems within 
this wider context. Further, we have analyzed several other quantities absent in Rcf. 
17. See, for example, Section 4.1. Of course, where we do have overlap, the results 
of 17 are reproduced exactly. 

2. Formulation of the problem and construction of the probability dis- 
tribution 

Consider a d-dirxiensional hypercubic, infinite lattice. Assume, that every lattice 
site is occupied by a single "white" particle (or ball) except for two sites: one of 




Fig. 1. The infinite sea of particles with one tagged (the solid circle) particle, and one empty site. 

them is empty (called "hole" or vacancy) and the other is occupied by a "black" 
colored ball, or particle, see Fig. 1. 

The system evolves according to the following rules: 

1) There is only an excluded volume constraint on the balls; no other interactions 
are present. 

2) The balls can take one step (equal to one lattice constant) jumps and only onto 
the empty site. 

3) The probability that the empty site will be occupied by one of its nearest neigh- 
bours is the same for every neighbour: p = \/2d. 

4) The black ball is distinguished from the others only by its color. Otherwise it 
obeys exactly the same rules l)-3) as the others. 

It is obvious that the hole performs random walk of Polya type BEm. Due to 
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rule 2) , the black ball performs a "passive" walk, moving only as a result of being 
"kicked" by the vacancy. One might as well say that it is driven by the Brownian 
vacancy. Our main concern is the motion of the colored ball. More precisely various 
properties of the walk performed by the black ball will be derived rigorously. 

In the following we write down the full probability distribution for the pair: 
vacancy and tagged particle. We have also found an alternative, and mathemati- 
cally more rigurous derivationa than the one presented below, however the version 
presented here is much shorter, and perhaps more intuitive, too. 

Without loss of generality let us assume that, initially, the black ball is in the 
origin, and the hole is at tq (rg ^ 0) (rg £ Z d ). After n steps by the hole, the 
black ball is to be found at site p and the hole, at site p + r, see Fig. 2. (p , r 6 Z d , 
r ^ 0). For later convenience, let us denote: R = p + r. 




Fig. 2. Initially the tagged particle (the solid circle) is in the origin, and the vacancy is at site rn 
(the empty square). After n steps the tagged particle is at site p and the random walker at site 
R = p + r. 

Let us introduce the following notations: 

• P n (s | sq) is the probability that the vacancy is found at site s after n steps, 
given it started at sq, n > 1, and s,sq G Z d . For n — by convention, 
P (s|s ) = S s ,s . 

• -P,t(s I sq) is the probability that the vacancy is found at site s (s,sq ^ 0) 
after n steps in the presence of an irreversible trap at site 0, given the walk 
commenced at site sq. In other words, P^(s | sq) is just the probability to get 
to s from Sq in n steps without having visited the origin. In the literature, 
the P^(s | sq)'s are called 'taboo' probabilities, see Ref. 6. Consistently, 
P t (s'|s) = fe,s. 

• 7r^(p, r 1 0, tq) is the probability that the black ball was driven k steps from 
to p, during which the random walker made n steps from tq to r + p. 

(s o ) 

Furthermore, will mean summing over all the nearest neighbours of site sq, 
s 

and in ^2 the summation is taken over all the lattice sites, 
s 
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Let us consider the referrence frame connected to the tagged particle. In this 
frame the walk performed by the vacancy away from the origin is just the regular 
Polya walk, but it can never visit the origin. Instead it "jumps over" the origin, by 
moving two lattice constants. This jump corresponds to the simple exchange of the 
vacancy and tagged particle in the original frame. With this observation, 

(0) (0) (0) oo oo oo 

n k n (p,r\0,r )=p k J2J2-Yl E E - E ^ElU^'EV' 

v\ V2 v k no=l«l=l nk=0 

Pi ir\-v k )Pl^_ x {u k | -v k ^)...Pl_ x {vs | -vjPl^ |r ) , (1) 

k>l, r o ,r^0. 

We had to sum over all the partitions of the time n into k subsets, as well a 
sum over all the possible directions (in number of 2d) from which the black ball can 
be "hit" . Obviously if i/j is a nearest neighbour of the origin, then — i/j is another 
nearest neighbour on the other side of the origin, straight accross. 

The k = case has to be treated separately: 

7r«(p,r|0,r )=^ pt(r|r ). (2) 

By analyzing ([!]), we infer: it k (p, r | 0, rg) = 0, if k > n, as it should be, since the 
black ball cannot be hit more times than the number of steps taken by the random 
walker. 



3. Generating functions 

Usually, it is more convenient to work with generating functions, instead of the 
distributions, themselves. In general, if the generating function for the distribution 
G n , i.e., G(£) = X)^Lo£™^™ ^ s known, the G„'s are retrieved by the complex 
integral along the contour C encircling the origin (counterclockwise):: 

c 

One can also quote the following relations, to be used extensively in our treatment: 
X> = ^/ 4^), and f> = G(l). (4) 

3=0 J c * ^ 3=0 

In the case of the Polya walk, the generating function is expressed as B0: 

where the integral is taken over the first Brillouin zone [— 7r, 7r] , and ui is the struc- 
ture function of the walk w(l) = ^ Ej=i cos (f?)> — n — h ^ w (1 = J2j=i h e j with 
{ej} as unit vectors in M d ). 
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For the simplicity of writing, let us introduce the following notations: 

t = t(0=P(0\0;0, (6) 

u = u(0 = P(a\0;0 , (7) 

h = h(t)=P(a\-ar,€), (8) 

b = &(£) = P{(3 1 -a; £) , /3 ^ ±a , (only for d > 2) (9) 

where a and /3 are nearest neighbours of the origin (| a |= 1, |/3|=1). t is the 
generating function for the return probabilities, u is for the probabilities of return 
to a nearest neighbour, etc., see Fig. 3. 




Fig. 3. The generating functions t, h, b and u. 

The connection between the taboo and Polya walk probabilities is most simply 
given in terms of generating functions!: 

P t (s|s ;0-P(s|s ;0-^(s|0;0 J P(0|s ;C), s,s ^0. (10) 
Next, consider the identities^!: 

^e ils = (2n) d S(l) , s eZ d , le[-7r,^] d , (11) 
s 

(0) 



Ve^=- W (l), le[-7T,n] d , (12) 
P 

7T 



In view of the above, the following useful relations valid for any 1 g [— 7r,7r] d and 
s^, sq 6 Z d , are easily proved : 

X> ils P( S |s ;0 = (U) 
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„ K J (27r) d 1 - &{X) 



(15) 



By inserting 1 = in ([T^) and using (|s|) and Jl3|), is obtained: 



(Si 



E P(a | s ; = ^ [P(s x | Sq; f ) - 5 Sl , Sq 



(16) 



A further consequence is with = (a is a nearest neighbour of the origin) and 
sq = in (|l^). Then, in virtue of notations (|^-0) the left hand side becomes 2du, 
whereas the right hand side (t — l)/p£- Thus: 



t - 1 



(17) 



Consider now = and sq = (3 be a nearest neighbour of the origin. Then from 
©: 

(0) 

^2P(a\0;€) = t + 2(d-l)b + h = ±. (18) 

a 

Let us now assume the notations: will mean summing over all the sites 

s 

(S]_) 

except for s = 0; ^ ' will mean summing over all the nearest neighbour sites of sj 
s 

except for s = (if that is the case). 

Similarly to Eqs. (|I^-|l5|) one can write : for any 1 e [—• 7T,7r] d and s^,sq e Z d , 
s ^0: 



EV ls pt (s | So;0 = _ 



e ils o 



Mi) 



ils 

1-— P(0|s ;0 



(19) 



i - 



-^(o|s ;0 



19J) and (g0|) can easily be proved using (|10|) and (y),(|15j). 
A list of useful particular cases: 



5> lls pt (s | a;0 = _ 
s 

E'pt(s|s ;0 = T 1 



-ila 



1 - V 1 " 



£wfl) I t 



1 - 



^(o|s ;0 



(20) 

(21) 

(22) 
(23) 
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E^lso;^^^, (24) 



and 

(0) 

£pt(/3|a;£) = -~, (25) 

with a a nearest neighbour of the origin and sq ^ 0. 

Let us now return to our original problem, and trade the discrete time variable 
n for the continuous, complex variable £, via: 

oo 

n fe (p,r|0,r ;0 = £ fV^p, r 1 0, r ) . (26) 

From d) and (|) and using the generating function for the taboo probabilities 
one finds: 
k > 1 : 

(0) (0) (0) 

n fe (p,r|0,r ;0 = (K) fe EE-E 5 * P f (r | - £)^ f (^ | -i*-* ; 0- 

...Pt^l-i/^^pt^lro^) , (27) 

jfe = : 

n (p,r|0,r ;0='5p,o J P t (r|r ;0- (28) 

In Appendix A it is shown that the distribution 7r^(p, r | 0, tq) is normalized. 
(This also serves as a check for the correctness of ([!]) and (g).) There we first show 

oo 

E E E%P+r^ (P, r 1 0, r ) = P„ (R I r ) . (29) 
fc=o p r 

so that normalization immediately follows Er^^I 1 "!]) = -0- 

In many cases, the integral in (||) cannot be carried out explicitly. Only asymp- 
totic behavior can be given in some limit. Of course, the limit of most interest here 
is 7i — ► oo. For this_purpose, we quote an extremely useful theorem, the Discrete 
Tauberian Theorem 10: 

The Discrete Tauberian Theorem: The relations 

oo - 

J2 GnC ~ (1 - O^— -) as £->l" (30) 

n=0 ^ 



and 



G n ~ —. . L(n) as n — > oo (31) 
T(k) 
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are equivalent, provided that n > 0; G n > (for large enough n); the sequence {G n } 
is monotonic (for large enough n); and the function L is slowly varying in the sense 
that — > 1 as x — ► oo for each fixed, positive A. 

4. Properties of the lattice walk performed by the tagged particle. 



4.1. The probability that the tagged particle suffers k encounters during 
n steps of the vacancy 

Let us denote by 7T„(rr)) the probability that the tagged particle has received k 
hits during n steps of the vacancy, given that the vacancy commenced at site tq 
(rn ^ 0). One can certainly write: 



^n( r 0) =EE' f »(''' r l ' r 0) ' 



(32) 



so that the corresponding generating function li fe (rQ; £) = Y^=o ^ri( r o)£ n becomes 

n^r ;0 = ££'n*(p,r|0,r ;0. (33) 
p r 



Inserting the expressions fl27|),(|28D into ( |33| ) and then employing (|23j),(g2|) and (|25|), 
we can perform the sums over r and the Vj-s. The results are 
k > 1 : 



: 



n fe (r ;£) = — 



n°(r ;0 



u 

1 

t 



fe-i 



^(0|r ;e) 



1 



1 



^(0|r ;e) 



(34) 



(35) 



Eq. ([35]) is the generating function for the probability 7r°(rg), that the vacancy 
avoids the origin during the first n steps. For k > 1, the explicit evaluation of 
the inverse transform (^ is possible only in one dimension (see Appendix B for 
the expressions of the corresponding generating functions). However, we do not 
find these results particularly illuminating and will not present it here. Instead, we 
analyze these probabilities in the n — > 00 limit, giving the leading behaviour in n 
through The Discrete Tauberian Theorem. 

d = 1 Using the expressions for the generating functions in d = 1 from Ap- 
pendix B and applying the Discrete Tauberian Theorem for k > 1 with k = 1/2, 
and for k — with k = 1, we find 



^(r ) 



1 



2 e -V2- 

IT 



l 



7^ and for k = : 7T°(r ) - 1 - e (36) 



l 0' 



d = 2 For simplicity, we shall treat here the case tq = a, only, i.e., the vacancy 
starts from a nearest neighbour of the tagged particle which is in the origin. As we 
will see, the probability for the tracer to suffer k hits scales as Inn. Thus we may 
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expect that, apart from the initial "wait" until the walker reaches the tracer for the 
first time (i.e., if tq is not necessarily a nearest neighbour), the leading behaviour 
to be independent of which direction the vacancy first hits the tagged particle. 



Thus, we take P(0|rQ;£) = u. Considering the limiting behaviour (B.10) for t 
and applying the Tauberian Theorem with k = 1, we arrive at: 

7r£(a) ~ r^- e"^, k > , n»l. (37) 
Inn 

this result shows that the number of hits the tagged particle received scales like the 
Inn, with an exponential (or geometric) scaling function. Note that, from Section 
4.2, distances for the Brownian vacancy driven walk do scale as the square root of 
k, the number of steps taken by the tagged particM Of course, the scaling function 
there is not Gaussian. 

d > 3 t(l) is no longer divergent, but is a finite number larger than 1. One can 
easily see that for 'n = oo': 



^ (r o ) = t(T) 



i 



P(Q\r ;i) 

t(l)-l • 



1 



1 m 



(38) 



If 



If rg is a nearest neighbour of the origin, then ^^(a) = 
R(0) denotes the probability of return to the origin, then according to Ref. 6 
R(0) = 1 - l/t(l), thus ft*, (a) = [1 - R{0)]R(0) k . Observe that, independently of 
dimensionality, the large fc-behaviour is a decaying exponential. 

Another interesting question is: 'what is the average number of hits suffered 
by the tagged particle during n-steps of the walker i.e., (fc)„(ro)'? Naively, one 
expects that it is simply the average number of visits the vacancy pays to the 
origin, when started from site rn: (7) n (rn). However, on closer examination, we 
see that the tagged particle wanders away from its original position (the origin) 
and towards the initial position of the walker, so that it is not quite at the origin 
all the time. Therefore, (/c) n (rg) can differ from (7) n (rr)). We shall compute the 
difference (fc) n ( r o) — (-0™( r (j)- ^ ms ^ s easu Y given as: 

(*> w (ro) - (J)„(r ) = ±-J ^ l-\l- t P (°l r 0^) • ^ 

r 

In one dimension, the large n-behaviour turns out to be: 

|r l+i 



(40) 



i.e., in the n — > oo limit the tagged particle receives 0.5 more visits than the origin. 

In two dimensions, after evaluating the cut in the integral ( |39| ) in leading order, 
the following behaviour is found: 

, > , , 21 r/x , 2 Inn , . 

(k „ r - 7 n (r ~ — - V (n ) + c] ~ — , n» 1 , 41 
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where ip(n) is the Euler psi- function and c is the Euler constant. After taking the 
limit n — > oo, one has (k) 00(1*0) = COoo(ro)- 

In d > 3, there are no singularities in the integrand and the difference approaches 
to zero rapidly with n. 

Observe that (fc)„(ro) > (J) n , i.e., the tagged particle receives more visits than 
the origin in any dimension. This is explained by the following. As it will be shown 
in Section 4.3, the center of the walk performed by the tagged particle is always 
shifted from origin towards rn (by an amount (p)n( r n))- Since rn is the center 
of the walk performed by the vacancy, it means that the average position of the 
tagged particle is always in between the average position of the vacancy and the 
origin. Thus, the tagged particle receives more visits than the origin. 

4.2. The probability that the tagged particle arrives at site p on the n-th 
step of the vacancy 

Let us denote by 7r„(p | rg) the probability that the tagged particle arrives at 
site p on the n-th step of the vacancy. Clearly, we have 

00 

^(p|ro)=EE V «(^ r l°' r o)- ( 42 ) 
k=o r 

while the corresponding generating function is 

00 

n( P 1 r ; = E E' 11 " to r I °> r o ; • ( 43 ) 

fe=o r 

In this sum, the k = term has to be treated separately, and it is: 



l-ip(O|r o ;0 



(44) 



Let us denote by ITi(p|rg;£) the rest (k > 1) of the terms in (|4^). Using (£3|), the 
sum over r can readily be performed to give: 

_ , oo (0) (0) 

ni(p|r ;0 = T ^E(^) fc E-E 5 * P ] Wu\-u k ^;£,)...P\u x \v ;0- 

3=1 

(45) 

After introducing the function 
(0) (0) 

- E^)" E - E e- 1{Ux) " k PH»k I -u^^e-^-K.e-^P^^ \ r ; 

fe=l V\ Vk 

(46) 

and using the integral representation of the Kronecker delta, one finds: 

ni(p|r ;0 = 1 r ^ / ^e-SW- (47) 
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Let us now evaluate the function 'I'j(x), regarding 1 as a parameter. Define first 
the following integral operator: 



(2n) d 1 - 



The result of the action of fj on a function of variable i? is a function of a 
variable x an d parameter 1. The k = 1 term in ( f46| ) is 

fJ x) (tfL4(r ;£,tf), (49) 

where 

^(r ;e^)=e^ r O-ip(0|r ;e). (50) 

After performing the sum over in the rest of the terms (k > 2), and comparing 
the resultant sums with (EH), we find the following integral equation: 

- f j x) (^M(r ; 6 #0 - j*l(0)f| x) + fW^)*!^) , (51) 

where Ij is just fj {^-l) applied on the constant (= 1) function. Substituting 
the right hand side of ( J5l| ) into the ^(i!^) on the right repeatedly, we arrive at: 

*!(*) = fW^jf { #l) (^(r ;^<?) - ^(0)1^ ^Jrf l} , (52) 

where 

f^^o^i^c^) + fj x) (^) + fj x) ( 1 ? 2 )f} , ' 2) (i? 1 ) 

+ fW^jf^^jfW^)..., (53) 



with 



= / ^^-X) - (54) 



Choosing x = in (153): 



i + fT 1 (U) (t9 1 )r}' 5l) 



In the following we investigate the operators T and T. By virtue of (12) and (|4q), 
the action of f on a function / : _ZR d — > (F may be written as: 

P — 7T P 
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with /ig as 2d, 1-dependent complex numbers. 

Since the /3's are lattice sites (elements of Z d ), they are of the form: 



/3=(0,...0, a ,0,...,0), a €{-1,1}, 

/i€{l,2,...,t0 



T 



(57) 



(/i labels the axes). Then, instead of (3 one can use the integer j given by the 
following one-to-one correspondence. 



3 = (1- o-)- +/i , 



(58) 



if /3 is known. Conversely, 



(7=1 — 2 



, fi=j-d 



(59) 



when j is given. An immediate consequence is that when c = +1, we have j € 
{ 1, d} and when <r = — 1 , we have j € {d + 1, 2d} . 

Let us denote by V(M d ) the set of functions v : JR d — > (Z* spanned by the basis: 



(60) 



where /3(j) is the nearest neighbor of the origin corresponding to the integer j via 



( p9[ ) and (57). (V(M ) forms a linear space with the usual rules of addition and 
multiplication of functions over the division ring of the complex numbers.) Note 
from ([36]) that the result of the action of f on any function / : M d — > (F is a function 
from V(M d ). Furthermore, since the operators T and T are linear operators they 
admit a representation by 2d x 2d matrices in this linear space with the basis as 

2d 

in (|60|). Let Tj and Tj denote the two matrices. Every element v = is 

3=1 

represented as a 2c? x 1 columnar matrix [w] in this basis: 



V2d 



(61) 



Let us now define the matrix P = [pj 1 ] , j, m £ {1, 2d} in this linear space, with 
elements given by: 



P " l = P(f3(j)\- a (m);0 



(62) 
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Explicitly: 



h 


b . 


. b 


t 


b . 


. b 


b 


h . 


. b 


b 


t . 


. b 


b 


b . 


. h 


b 


b . 


t 


t 


b . 


. b 


h 


b . 


. b 


b 


t 


. b 


b 


h . 


. b 


b 


b . 


. t 


b 


b . 


. h 



(63) 



Next, let us consider the diagonal matrix Lj 



L 



, with matrix elements: 



jm _ e -ilf3(j)fim 



(64) 



where S" 1 (Kronecker delta) are the elements of the 2d x 2d unit matrix 1. 



The operator f will be represented by the matrix Tj 
d d $ e i^M+gO)] _ a _ 



, with elements: 



From the above, we have: 



= KP(/3(j)|-a(m);0e 



(65) 



T X = K^P . (66) 
Then, based on (E^) , it is easy to write down the representation for the operator T, 



tV 1 
L h 



t { = i + r t + rf + rf + ... = (1 r,) 1 



(67) 



Examining (55), it is obvious that we need to calculate the action of the operator T 
on two functions: the constant (= 1) function and e 1 r 0. If [q] denotes the vector 
representing the action of T on the function unity and [rj] represents the action of 

f on e 1 r we have: 



[4] = Vp^uLi^i] , 
[t 1 ]=pZL 1 [Yo], 



(68) 
(69) 



where [4>\\ is the eigenvector of P corresponding to the eigenvalue Ai = uj (pt;) (see 
Appendix C) and [Yq] is the vector: 



[Yo] 



P(/3(l)|r ;0 
P(/3(2d)|r ;0 



(70) 
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Note that 



Since 



[faf[Yo] 



P(0\r ;0 



2d 



3=1 
2d 



3=1 



and V (0) — 1, we arrive at: 

vWxpT^yo] - £?p(o I r ; O[0i] T T 1 L 1 [0 1 ] 



*l(0) = 



1 + 



(71) 

(72) 
(73) 

(74) 



In Appendix C, we show how to diagonalise the matrix P and write it in the form 
P = C _1 r2C, with $7 a diagonal matrix containing the eigenvalues of P on the 
diagonal. The columns of the matrix C _1 are the right eigenvectors of P. Therefore 



i 1 ] T C- 1 = [10 ... 0] , and C[fa] = 



(75) 



Thus 



[fafT^fa] = [0 1 ] T (l-KLiP)- 1 L 1 [0 1 ] = [^HL^-KP)- 1 ^] = 
[0 1 ] T C- 1 (CL_ 1 C- 1 -ptn^Cfa] , 

where we have used the fact Lj -1 = L i. Therefore [0i] t TjLj [fa] is nothing but 
the m = 1, j = 1 element of the matrix (CL_jC^ 1 — p^fi) -1 . 

For I = 0, we have Lq = 1, thus the first element in the above matrix is: 



[^] T T o Lo[0i] 



1 



1 



1 — p£,\i 1 — u 

With the help of ([FT]), the first term in the numerator of (|7l| ) becomes: 

1 1 



Vp£[</>i] T ToL [y ] 



1 -u t 



and the second: 



u 1 



f- P(0 1 r ; [fa] 1 T L [fa] = £- — P(0 1 r Q ; £) 



Collecting all the above together, we have: 



M0) = T ^u]p(O\v O ;O 



(76) 



(77) 



(78) 



(79) 
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One check for the above calculations is to see if the distribution 
ized. For this purpose, we consider: 



is normal- 



^n(p|r ;0 = T 1 - 



l-ip(0|r ;f) 



X) n i(pl r o;0 



1 



l-ip(0|r ;0 



1 - u/f 



*o(0) = 



(80) 



showing that normalization is indeed properly satisfied. 

Considering the first term ( Ejj ) and the expression (fn^) , the generating function 

is: 

7T 

n(p|r ;0 = T ^/ 7 ^ e - ilp {l+(l-^)[*l(0)-* (0)]} , (81) 



where 'I'j(O) is given by (|74|). 

In the following we examine the meaning of formulas (|8l]) and (|74| ) , by consid- 
ering them in different spatial dimensions. 

d = 1 . This case is considered more as a check for the validity of our formulas since 
in one dimensions the problem is rather trivial: the random walker is either on the 
left or on the right of the tagged particle, and the only way to change sides is to 
exchange positions with the tagged particle. So, if initially the tagged particle was 
in the origin and the walker was to the right (tq > 0), the only two sites visited 
by the tagged particle are sites and 1. Nevertheless, one can ask a nontrivial 
question: 'what is the probability to find the tagged particle on sites and 1 after 
7i-steps of the random walker, given the walker started from site ro ?' This can be 
answered in the framework of our theory. Without loss of generality, we assume 
r > 0. 

The matrix P is in the simplest possible form: 



h t 
t h 

The computation of 'I'i(O) is rather straightforward using ( |74] ) 

t l-e-"\ 



*j(0) = * (0) 1 - 



1 + 7 



where *o(0) is given by (ff9|) and 

P(0|r ;£)=«(~) 



i-o-l 



u 
1 



See Appendix B for the generating functions. The integral over I in (p 
and we finally obtain: 



(82) 

(83) 

(84) 
is trivial 
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Indeed, as expected, there are only two sites the tagged particle may visit! Moreover, 
using Eqs. (||) and (||), we have Woo(p\fo) = (l/2)(^,o + <5p,i) ■ I n other words, the 
tagged particle spends half of its time on site 1 and half on site 0. The n-dependcncc 
is found by simply performing the integral with respect to £ . 

d = 2 . This is the nontrivial case. In the theory of recurrent events it is shown 
under very general circumstances that the return probability R(sq) to a site sq 
for a lattice walk, is simply expressed by: -R(sq) = 1 — 1/P(sq | Sq; 1~), where 
P(sq I sq;£) is the generating function for the lattice walk (which can be almost 
anything, not necessarily Brownian) . In Ref. 9 we proved that the same conditions 
are applicable to the walk performed by the tagged particle as well, provided we 
replace P„(s | sq) by W n (p | rg). Here, n means both 'time' and the number of steps 
taken by the vacancy, p "corresponds" to s and, according to our initial conditions, 
sq is replaced by while rg serves as a parameter. Therefore, the walk performed 
by the tagged particle is recurrent if II (o | rg; 1~) is divergent, which is indeed the 
case, see below. 

In two dimensions our matrices are four by four, and the process of evaluating 
II(p | rg;£) via ( |8l] ) and (frij) is not particularly illuminating. We report only its 
final, analytic expression: 



n(p|r ;0 



P(0|r ;£)] +I -— £>(p;0£i(ro;0> 



(86) 



where 



3j(p;0 



(2tt)2 



25(1) 



3 = 1, 



(87) 



and 



/Ci(r ;0 
^ 2 (r ;e) 

£s(ro;0 
^ 4 (r ;0 



^P(0|r ;0, 

|[p(/3(l)|r ;e)+P(-/3(l)|rg;0 
-P(/3(2)|r ;e)-P(-/3(2)|r ;£)" , 

P(/3(2)|r ;£)-P(-/3(2)|r ;e) 

P(/3(l)|r ;£)-P(-/3(l)|r ;£) 



2V2 



(89) 
(90) 

(91) 



The quantities (1) are proportional to the elements (up to a common factor) of 
the first row of the matrix (CL_jC _1 — p^Q)^ 1 and they are all expressed in terms 
of V = V(0 = |[|E(0 - 1 - (1 K(0] and W = W(£) = |[|E(0 - 1] which 

are the second and third eigenvalues of matrix P, A2 and A3 up to a factor p£ = £/4, 
Here, K(£) and E(£) being the complete elliptic integrals of first and second kind, 
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respectively (see Appendices C and B). 



Ai(l) = -{1 + VW)W + (1 + VW) + {V + W)W - cos(h) + cos(l 2 ) 



-(V + W)cos(h)cos(l 2 ) 



A 2 (l) 
A 3 (l) 
A 4 (l) 



{l-W 2 )-\cos{h) + cos{l 2 ) 



-=sin(l 2 ) (1 + VW) - (V + W)cos{h) 
V2 L 

— l T sin(li) \{l + VW) - (V + W)cos{l 2 ) 
V2 L 



(92) 
(93) 

(94) 
(95) 



and 



D(l) = 



(1 + VW) (l + W^j [l - (A + B)~[cos(h) + cos(l 2 ) 

+ABcos(h)cos(l 2 )} , 



with 



A 



W + j 
1 + W± 



D 



V + W 
1 + VW 



(96) 



(97) 



Formulas (|8q-|97|) give the exact expression of the generating function n(p | tq; £). 
From (^6|) we see that the functions gj carry the information on the spatial compo- 
nent of the walk performed by the tagged particle while the functions JCj describe 
the dependence on the initial condition of the vacancy. Rewriting the cosines and 
sines in formulas ( |9^ - 195| ) in terms of exponentials, it becomes clear that one needs 
in general to evaluate the £ — > 1~ behaviour of integrals of type: 



r(s;0 



d 2 l 



(27T) 2 D(l) 



(98) 



By further introducing the vectors 



[T(p;Ok 



r(p + /3(l);0 
r(p + /3(2);0 
r(p-/3(l);0 
r(p-/3(2);0 



r(p + /3(l)+/3(2);£) 
r(p-/3(l) + /3(2);£) 
r(p-/3(l)-/3(2);£) 
r(p + /3(l)-/3(2);C) 



(99) 



characterizing the nearest neighbour and next to nearest neighbour terms corre- 
sponding to site p, the Qj functions are simply expressed as: 



+ VW)w(l-j)r(jr,£) + 

(i + vw) + (v + w)w] (i - 1) W T [r(p;0] 7 



-(v + w) i-~ )[<t>iY[r(p;0] 



(100) 
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9z(p;0 = 


- 


-W 2 )(l- 


l) 


fc] T [T(p;0]„. , 


9a(p;0 = 


- & 


-VW) (l - 




[^3] T [r(p;0]n.+ 






f W) (l- 


?)( 


'W T + [^] T )[r(p;0]n.n 






-VW) (l - 


-?) 


[04] T [r(p;O]n. + 




+ \ {v - 


f W) (\ - 




W T )[r(p;0]n.n 



(101) 



(102) 



(103) 

The £ — > 1 _ behaviour is determined by the singularities of -D(l). Analyzing the 
functions ^4 and B one can rigurously show that there is only one zero of D(l) in 
the first Brillouin zone, namely 1 = 0. 

For s = 0, can be evaluated exactly: 



m = I ^TvwWTwf) T~ab K ( 1^) • « 

For s / the exact evaluation of the integral is cumbersome (it can be given 
recursively and in terms of Legendre functions of second kind) and it is not our 
aim here. Instead we analyze the £ — * l - behaviour. In this limit u/t — ► 1, 
V ->V = ±-1 ~ 0.273, W -> W = 1-1 ~ -0.363, A^A = 1, £ -^B ~ -0.100, 
see also Appendices B and C. The limiting behaviours for both s = and s ^ 
cases are computed to be: 

m o - £ (7r _4-i) K ( x - ^ - «/*)) - ^ 1_ • ( 105 ) 

and 



F(S; C) ~ T (7T - 2)(7T - 1) ^ i 2 ^ 31 I S I ^ 1 ' ^ 1_ > S ^ - 

(106) 

with Kq being a modified Bessel function. Let us now keep | s | constant. Since : 
K (z) ~ — lnz as z — > and K(z) ~ | In ( tzt) as z — > 1~ it follows that both T's 
diverge and have the same leading term, namely: 

r ( s; ^i ( ,-4-i) h (R)' (107) 



Taking into account that, in ( pM[p^ ), [c^-] T [r(p; £)]„. and [</>,] T [r(p; £)]„.„. for 
j > 2 involve differences of T's, we see that the nearest neighbour terms and the 
next nearest neighbour terms (and therefore <72, 53, gi) are negligible compared to 
the terms in gy. Thus, to leading order, we have: 

fttat )-(i-!)(i-SW i4t)+o(i-t) ■ f-i". 
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ffj(p;0 = + O(l- j) , e-l", j = 2,3,4. (109) 

One can see that the leading one is gi(p;£). Although the generating functions 
P(s I rg;l~) are divergent in two dimensions, the quantities /Cy(rg;l _ ) are all 
finite numbers. This is the consequence of formulas ( p8[]9l| ) and the fact that when 
£ — > l - the probability generating functions for the Polya walk all diverge in the 
same way, more specifically lim [£(£) — P(s | rg;£)] is a finite number, for any s 

and tq (see for e.g. Ref. 6, p. 145). Thus /Ci(rg; l - ) = 1. Using this, it also can be 
proved that the k = term ([II]) is negligible, compared to the giK,\ term, in this 
limit. 

Collecting the above, our conclusion is: when p is kept finite, the generating 
function Il(p|rQ;£) is, to leading order, 

n(p|r ;£)~ (l-z) j^—Jl -7)111 (j^—), £-1-, P finite. (110) 



tt / 1 - £ \ t) V 1 - t 
Using the Discrete Tauberian Theorem and the fact that 

u TT 

1 t r- , as < 1 - £ < 1 

1 In I -± 



we find: 



i-C 



_ 1 1 \ 1 .In (Inn) 

^n{p\TQ) ~ (7r - 1) — , 71 -> 00 , p finite. (Ill) 



From (11C), 11(0 | rg; 1~) is clearly divergent, proving our assertion that the walk 
in two dimensions is reccurent. The tagged particle is guaranteed to visit any site 
in the n — > 00 limit. 

Scaling. It is well known that, for the pure random (Polya) walk, distances scale 
as the square root of time, while the scaling function is, apart from a prefactor, 
Gaussian. However, for the Brownian vacancy driven walk, the scaling properties 
are quite different. When | p | and n are both large (keep tq constant), the only 
surviving term is again g\K,\. Since /Ci(rg; 1~) = 1, we infer 



e-fi-, ipi»i, (112) 



or 



7f„(p|r ) ~ 2(tt - 1)— K 2^tt{tt - 1)^= , n -» 00 , |p|» 1 . (113) 
Inn V ylnn/ 

This is the scaling function in the corresponding limit. As we see, instead of being 
Gaussian, it is the modified Bessel function Kq. From its argument, it follows that 
distances scale like square root of the log of time. 

In the following Section we calculate the first and second moments for the 
displacement of the tagged particle. However, instead of using the distribution 
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7r„(p|rg) derived in this Section, we will use a slightly different, but easier, way to 
arrive at the desired results. 

4.3. Average displacement and standard deviation of the tagged particle. 

Let us define the first and second moments associated with the tagged particle's 
displacement p up to time n: 

oo 

(p)„(r )^ /) f„(p|r )=^ /) ^^4(p,r|0,r ) , (114) 
p p k=o r 

and 

oo 

(P 2 )n(r ) =£AW(/»l r 0) =EP 2 EE' ^(P,r|0,r ) , (115) 
p p k=o r 

respectively. The first of these is, of course, the average displacement for the particle: 
(p) n (rg). The second is used in finding the "spread", i.e., the standard deviation: 



2? n (r ) = ^/<P 2 )„(r ) - [<p)„(r )P . (116) 
Instead of computing directly the sums above we turn to generating functions 

oo oo 

A(r ;0 = ^<p)„(r )r = EE^E' n fc (p,r|0,r ) , (117) 

n=0 fc=0 p r 

and 

oo oo 

S(r ;0 = E(P 2 )«( r 0)r = EE^E' n fc (p,r|0,r ) . (118) 

n=0 fe=0 p r 

By peforming the sums over r and p we arrive at the final expressions for A(rg; £) 
and 5(rg;£) (see Appendix D): 

A(ro;0 = Y^Trb M(t ° ;0 ' (119) 



1 — 4 1 — u/t 1 — q t 



(120) 



where 



and 



(0) 

M(r ;0= P eE^|ro;0, (121) 



q=rf(h-t). (122) 



Note, for rg = a in ( 121 ) we have: 

M(a;0 = -?«- (123) 
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The corresponding moments will then be expressed by: 

M(r o ;0 



i i di i mo 



(124) 



and 



r 



>™ = eTTi-e t^a t - g (125) 

r 

There is a simple expressionEl for M(rg;£) showing that it is in fact a vector 
pointing along rg: 

€ 

M(r ;0=r ||p(0|r ;CK- (126) 
o 

This means, that the average of the tagged walk is always shifted towards the initial 
position of the vacancy. 

In the following we extract the n — > oo behaviour for the different cases. 
d — 1 Without loss of generality, assume ro > 1. Then, from Appendix B: 



M(r ;0 = ^ e i . and 9 = -^^jE . (i 27 ) 

Therefore, for the average displacement, we find: 

]—) - (128) 

Using (||), we find immediately (p) 00(^0) = 1/2, which is certainly correct because 
half of the time the tagged particle is at the origin, and half of the time is on site 
1. The convergence to 1/2, i.e., the n-behaviour in leading order is given by: 

(P) n (r ) - \ (e^^) 6l . (129) 



Observe that, as pointed out below Eq. (|126j) on average, the tagged particle suffers 
a finite displacement toward the initial position of the vacancy. One can consider 
this as an effective attraction between the tagged particle and the random walker. 
This is property valid in any dimension. 

The second order moment turns out to be: (/9 2 ) n (r ) = (p)„(ro). Thus, the 
spread is Poo^o) = 1/2, as expected. 

d = 2 To obtain explicit expressions for the general case is non-trivial. For 
example, if the tagged particle and the vacancy are initially very far apart, we 
expect that the particle will experience the first hit to come from all directions 
with equal probability. Then its average displacement will vanish. So, we choose to 
present only the case where the vacancy-tracer is initially a nearest neighbour pair, 
i.e., rg = a, where a is a nearest neighbour of the origin. 
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The average displacement is then given by: 

d£ 1 (~q) 



I \ ( \ 1 I 1 " 



(130) 



and in the limit n — > oo (see Appendix B for the expression of q): 

1 7T — 2 

(p)oo(a) = ol = 0.266529a. 

2 7T — 1 



Applying The Discrete Tauberian Theorem, the second moment is found: 
(P 2 )n(a) 

this leads to the spread: 

m 1/2 

1 



1 , 1 

— ; 111 U : , fl»l. 

n(n- 1) i-km ' 

*> ' T7.7T 



(131) 



(132) 



V n {a) 



1 . 1 fw-2 
inn 



7r(7T — 1) 



4 Vtt- 1 



Inn , n > 1. (133) 



i.e. it diverges as Vmn. This is exactly the result we anticipated from the scaling 



given by ( 113| ). 

d > 3 Since there is no divergence in this case, the first and second moments 
are finite even for n — oo: 



<p)oo(r ° ) = T^T ' (P 2 )oo(r ) = ^|^(0|r ;l). 



The spread, 



^oo(r )~|i^P(0|r ;l) 



1-9(1) 



1/2 



(134) 



(135) 



is also just a finite number. 

If the vacancy again started from a nearest neighbour a of the tagged particle, 
then, we have, via (|l23|): 



Woo (" J 



g(l) 
1-9(1) 



a , and 00(a) 



1 



1-9(1) 



1 - q(l) 2 t(l) - 1 . (136) 



In particular, in three dimensions t(l) = 1.516386 and q(l) = —0.209.. (see Ap- 
pendix B), so that 



(p)oo(a) = 0.1728.. a , and V^a) = 0.5549. 



(137) 



5. Conclusions and summary 

We investigated the walk performed by the tagged particle in the presence of a 
single vacancy in a d-dimensional, infinite square lattice. The only "active" object 
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is the vacancy, while the tagged particle performs a "passive" walk, being shifted 
only when it exchanged positions with the Brownian vacancy. We constructed the 
complete probability distribution, which is quite complex, incorporating three spa- 
tial variables (rg, p and r) and two temporal ones (n and k). However, physically 
relevant informations can be extracted from different projections of the full distri- 
bution. An example is W n (p \ rg), the probability distribution for the tagged particle 
to be at site p on the n-th step of the vacancy. Another is 7r„(rn), the probability 
that it has made k steps during this time. Closed form expressions for the corre- 
sponding generating functions, in terms of the generating function of the Brownian 
walker alone, were obtained. Exploiting the powerful Discrete Tauberian Theorem, 
we find the limiting behaviours of the distributions themselves, as n — > oo. We 
showed that the Brownian vacancy driven walk is reccurent in two dimensions only. 
Other interesting results associated with a tagged particle in d — 2 are 

(i) the number of steps it takes scales as Inn (as n — > oo) 

(ii) the limiting displacement- distribution is not a Gaussian but rather the modi- 
fied Bessel function Kq 

(Hi) its standard deviation scales as y/]nn, in contrast to the usual y/n scaling 
Brownian walks 

By computing the first moment for the displacement-distribution of the tagged 
particle, we find that there is an effective attraction towards the vacancy which 
shifts the average position of the tagged particle from its initial point. In d=l, 
2, and 3, this shift is, respectively, 0.5, 0.266529.., 0.1728.. in units of the lattice 



spacing. For higher dimensions this is given by the first expression in (134). We 
have seen, that an immediate consequence of this attraction is the tagged particle 
being visited more frequently than its original site. However the difference vanishes 
as n — > oo, for d > 2. 

In d > 3, the standard deviation of the displacement-distribution does not di- 
verge as n — > oo. Instead, it converges to a small number, e.g., 0.5549 in d = 3. 
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Appendix A 
Normalization 

When proving ( ^9j ) it is better to switch to generating functions. (^9|) is equiv- 
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alent to: 



E E E' 5 R,P+r n" (P, r 1 0, r ; = P(R I r ; 



(A.l) 



fe=o p r 



Separate the k = term in (A.l): 



T ° = E E^R.P + r n ° (p, r 1 0, r ; ^ P f (R | r ; 
p r 



(0) 



1 



P(R|r ;£)--P(R|0;£)i 5 (0|r ;0 



(A.2) 



As for the rest (p = R — r) : 



(0) (0) 



^-EEE^Rp + rn fc (P,r|0,r ;0^ / TCT^E^E'-E 

fe=i p r J v 7 fc=i i^i ^ 



E'e^Crl-^jO 



e-^pt^i.^^^.^-ilnpt^i^.^) . (A . 3) 



Using Eq. ( BjJ ) with a = — i/j., the term in square bracket is computed, and: 

ilR oo (0) (0) (0) 



T 



d a l e 1 



(27r) d 1 - 

pt^l-^.^O.-.e-^pt^l-iy^Oe-^P^^Iro;^) . (A.4) 



According to the square bracket in ( A.4 ), T splits in two terms (T = 81+62)- Using 
notation (E6J) , the second term becomes: 



it j d d l e ilR 
tj (2ir) d 1 - ^(1) 



*l(0) 



(A.5) 



while the first: 



.ilR 



(0) (0) 



(2tt)<* 1 - £u>Q) ^ 



E^ fe EE- 



(0) 



^pt^i-^o 



-ip+^-xl-i/fc-ajO-e-^PUi/xIrojO • (A.6) 



We again need to separate the k = 1 term in #1: 



<M e 



1R 



(0) 



{2n) d 1 - 



K E^Kko;?) +(^>2 terms). (A.7) 
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The term in square brackets in (A. 6) is computed via (f25|), while the term in curly 
brackets in (A. 7) is expressed by (|24|). Combining these, we finally get for the first 
term: 



h = ip(R|0;£)P(0|r ;£) ' ? " 



d d l 



,ilR 



(0) (0) (0) 



t J (2n) d 1 - 



e il» k -i pt (lyfe _ i | _ Uk 2 . £) . ,. e -il^ P t ( Vi | r() . e) 



fc=2 



^1 ^2 ^fc-l 



I u r d a l e'^ 

= -P(R|O;0P(O|r o ;0 + jjp^O). 

— 7T 

Adding the first (#1) and the second term (#2), is obtained: 



T = ip(R|0;£)P(0|r ;£) 



(A., 



(A.9) 



(the terms containing the \t-s cancel). A dd (|A.9|) to ( |A.2| ) and after cancellation, 
(A.f ) is gained. Therefore normalization (A.l) holds. 



Appendix B 

Generating functions 



For dimension one, the explicit expression for the P(0 1 R; £) generating function 
is well known Q: 

\R\ 



P(O|P;0=t 



with: 



Thus, according to (|) (or (|l7|) ), and (|): 

„,2 



U = t 



h = 



t 



In two dimensions 



IE 



t = P(0|0;0 = -K(0 



(B.l) 
(B.2) 

(B.3) 

(B.4) 



where K(£) is the complete elliptic integral of first kind. Then, according to (P_7J): 



1 

u = — 



-K(0 - 1 



(B.5) 



Quoting Ref. 15, the probabilities P(0 1 R; £) with R = R(m, to), me Z (R(to 1; to 2 ) 
is a lattice vector with components m\ and to 2 ) can be expressed in terms of Leg- 
endre functions of second kind: 

P(O|R(to,to);0 = ^Q m -l/2 ^ " 2 



e/2 



(B.6) 
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Note, that for the Polya walk P(0 | R(mi, m 2 ); = P(0 | R(mi, -m 2 ); £) = P(0 | 
R(— mi, m 2 ); £) = P(0 1 R(— mi, — m 2 ); £) for any mi, m 2 e Z. 



For m = we regain (B.4). For m = 1 
ft 



P(0|i2(l,l);0 = -^ 



1 - ^ ) K(0 - E(0 



(B.7) 



with E(£) as the complete elliptic integral of second kind. 
From (ELI) follows: 



2 4 
/i = -K(0 + -s 



-E(0 - 1 

7T 



Therefore, from (122): 



lE(f) - 1 



(B. 



(B-9) 



When £ — > 1 , t, u, b and /i all diverge, since K(£) has a logarithmic singularity in 
£ = 1. Accordingly, we have the following asimptotic behaviour @: 



7T 1 — S 



as £ 



r 



(B.10) 



Note that for $ = 1, q(l) = -(1 - 2/tt) = -0.36338. 

In dimensions three the value for t(l) is given analitically Q as t(l) = g^lX^j) 
r (^) r (^) r (|i) = 1-516386. We numerically estimated the value of 6(1) as 6(1) = 
0.331. Thus, from equation @ t(l) + 46(1) + h(l) = 6u(l), and one finds /i(l) = 
0.2579. Therefore q{\) = -0.209. 

Appendix C 
Diagonalization of P 

P is easily diagonalised. It has three different eigenvalues (for d = 1 only two, 
Ai and A 3 ): 



A 9 — A-i 



Ai 


= A a 


-x d 


= A 2 


A 2c / 


= A 3 



h + 2(d- 1)6 + * 
ft, - 26 + 1 , 



The corresponding system of orthonormal eigenvectors is: 

1 



1 



1 



1 



as the single eigenvector corresponding to Ai , 
1 



[0iF 



VW - 1) 



1 .. (J - 1) 



1 



1 



-1 .. (j - 1) 
3 e {2,3,...,d} 



1 .. 



(C.l) 
(C.2) 
(C.3) 



(C.4) 

- 1 
(C.5) 
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as the d — 1 eigenvectors corresponding to A2 (the small vertical line is between 
elements d and d + 1). The elements —1 in [<f>j] T are in the columns 1, d — j + 3, 
d — j + 4, d — 1, d and d + 1, 2d — j + 3, 2d — j + 4, 2d — 1, 2d, respectively, 
while the elements j — 1 are at the columns d — j + 2 and 2d — j + 2, respectively. 
And, finally, 



1 
71 







-1 







■•• 1 ••• 

ie {d+l,d + 2,...,2d} 



(C.6) 



are the d eigenvectors corresponding A3, with nonzero elements at columns 2d+l — i 
and 3d + 1 — i. It is certainly valid: 



] T [4>j]=S id , i,j£ {1,2, 3,..., 2d} 



(C.7) 



d = 1. The eigenvalues with the corresponding eigenvectors are: 

and 



2 1 

X 1 = h + t=-u, [<j>x] = —= 

? V 2 

2m r, ! 1 

\ 2 = h- t = -— , W = - 



-1 



1 



with t and u given by Eqs. (B.2) and ( |B.3| ). While Ai diverges as £ 
convergent (u/t 1). 

d = 2. In this case the eigenvalues are given by: 



(C.8) 
1~, A 2 is 



4 4 

Ai = h + 2b + t=-u = — 



-K(0 - 1 

7T 



A 2 = ft - 26 + t = — 
4 



c 2 



ft, - < 



i E (o-i-(i-a-K(o 

7T 7T 

"2 



E(0 - 1 



-W 



(C.9) 
(CIO) 
(C.ll) 



and the corresponding eigenvectors by: 





" 1 " 




' -1 " 


1 


1 




1 


2 


1 




-1 




1 




1 



1 

71 



= 



V2 



(C.12) 
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First and second moments 
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Using the expressions (|27|) and (|28|) for the generating functions, and performing 
the sums over r and p, we have: 



A(r ;0 = i^E(^) fc E-"E^ + "-+^) pt ^|-^--^"- pt (^|r ;e). 

(D.l) 



(0) (0) 



fc=i 



Using ( fL0| ) and the notation (121), the k = 1 term is 

(0) 



M(r ;0- 



(D.2) 



Defining: 



(0) (0) 

Ai(r ; = £(K) fe E -E^ + - + ^-O^K I -^^O-^K |ro;0 . 



fe=2 



and 



(0) (0) 

A 2 (r ;e)^E(P0 fc E-E^ pt KI"^--0- J P t (^|r ;0, 



(D.3) 
(D.4) 



fc=2 



(D.l) is rewritten as: 

A(r ;0 

It is easy to see that: 



1 - u/t 



[M(r ;0 + Aifo; + A 2 (r ; 0] • (D.5) 



Ai(r ;0 = 7T 1 -^7 A ( r o;0 



t 1 - u/t 

With the aid of formulas ( p5|) and ( |121| ) for A 2 (rQ;^) we obtain: 



A 2 (ro;0 = T A -M(r ;0 



(D.6) 



(D.7) 



Inserting (|DJ) and ( p/H ) into flDJ), (JTT^) is gained. 

Proceeding similarly as for the first moment, the generating function for the 
second moment becomes: 



%o; o = X -Y^r E(^) fe E - E(^+-+^) 2pt to I £)- pt to I r o;0- 



(0) (0) 



fc=l 



(D.. 



Separating the term fc = 1 in (D.8) and using the fact that (u 1 + ... + u^) 
(u ± + ... + + 2(u ± + ... + i/ fe _ 1 ) + 1 (u t 2 = 1, Vi), SXro;^) is rewritten a 



1 - u/t 



^(0|r ;0 



+ Si(r ;0 + 2S a (r ;0 + S3(r ;0 



, (D.9) 
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where 

(0) (0) 

^i(r ;0 = E(^) fc E-E(^ + - + ^-) 2pt KI-^-^)- pt (^l r o;0, 



k=2 v\ Vk 



(D.10) 
(0) (0) 

s 2 (r ; - X)(pO* E - E(^ + - + "fc-iK^fa I -v k -i;0...PH^ I r ; , 

fc = 2 Kl 

(D.ll) 

(0) (0) 

^(r ;0 = E(^) fc E-E pt K|-^-^)- pt (^l r 0;0, (D.12) 

fc=2 vi u k 

In virtue of p4[ ) and ( p5| ) is not hard to see that: 

Inserting these into (jDjjj), ( |l20| ) follows. 
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